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The spontaneous breaking of chiral symmetry implies the presence of a short-distance scale in 
the QCD vacuum, which phenomenologically may be associated with the "size" of the constituent 
quark, p ~ 0.3 fm. We discuss the role of this scale in the matrix elements of the twist-4 and 3 
quark-gluon operators determining the leading power (l/Q^-) corrections to the moments of the 
nucleon spin structure functions. We argue that the flavor-nonsinglet twist-4 matrix element, f2~'^, 
has a sizable negative value of the order , due to the presence of sea quarks with virtualities 
~ in the proton wave function. The twist-3 matrix element, d2, is not related to the scale p~^ . 
Our arguments support the results of previous calculations of the matrix elements in the instanton 
vacuum model. We show that this qualitative picture is in agreement with the phenomenological 
higher-twist correction extracted from an NLO QCD fit to the world data on g\ and p", which 
include recent data from the Jefferson Lab Hall A and COMPASS experiments. We comment on 
the implications of the short-distance scale p for quark-hadron duality and the a;-dependence of 
higher-twist contributions. 
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I. INTRODUCTION 

Polarized deep-inelastic scattering (DIS) has exten- 
sively been studied in fixed-target experiments with elec- 
tron and muon beams The main goal of such mea- 
surements is the extraction of the polarized parton densi- 
ties in the nucleon from the data in the region of approx- 
imate Bjorken scaling, where the Q^-dependence is gov- 
erned by perturbative QCD (DGLAP evolution). How- 
ever, most of the precise data are at momentum transfers 

~ few GeV^, where non-perturbative effects in the 
dependence (power corrections) generally cannot be ne- 
glected. Estimates of these corrections are needed in or- 
der to include these data in the QCD analysis. The effect 
of phenomenological power corrections on the extracted 
parton densities has been investigated in Refs. 

With the polarized DIS data becoming more and more 
precise, the study of the power corrections themselves has 
emerged as an interesting subject. The l/Q^-corrections 
to the lowest non-zero moments of the spin structure 
functions, gi and g2, are governed by matrix elements 
of QCD operators of twist 4 and 3, describing non- 
perturbative correlations of the quark and gluon fields 
in the nucleon ^|^: 

(iV|^7„gG'5V|iV) = 2hs^ (1) 

-I- traces, (2) 

where G"^ = {l/2)e°'P^^G^8 is the dual gluon field 
strength, and and s" the nucleon four-momentum 
and polarization vector; we follow the conventions of 
Ref. p 113. By extracting the coefficients of the l/Q^ 
corrections from the data one can thus obtain informa- 
tion about the structure of the nucleon in QCD. The 
SLAC E155X experiment and the Jefferson Lab Hall 
A experiment ISj have extracted the twist-3 matrix el- 



ements from combined measurements of the structure 
functions gi and 32, and found surprisingly small val- 
ues, d2 ^ 10~^. Recent analyses have also attempted to 
extract the twist-4 matrix element, /o, from the power 
corrections to the first moment of gi ^ 0, 0, ; see 
Ref. 113| for a review. 

The theoretical estimation of the higher-twist matrix 
elements Q and Q is a challenging problem, requiring 
a description of the nucleon in terms of QCD degrees 
of freedom. The key question is which non-perturbative 
scales govern the quark-gluon correlations measured by 
the twist-4 and 3 operators. This question is intimately 
related to the role of the "vacuum structure" of QCD in 
determining the structure of the nucleon. 

There is strong evidence for the existence of a short- 
distance scale in the QCD vacuum, significantly smaller 
than the size of the nucleon. It is determined by the char- 
acteristic size, p, of the non-perturbative field configura- 
tions instrumental in the spontaneous breaking of chiral 
symmetry. Numerous observations suggest a "two-scale 
picture" of hadron structure, 

p < i?, (3) 

where i? is a typical radius of the nucleon (say, the 
charge radius). Phenomenologically, the short-distance 
scale may be associated with the "size" of the constituent 
quark. In fact, the success of effective models based on 
constituent quark degrees of freedom could not be ex- 
plained without the hierarchy It is natural to ask 
what the existence of this short-distance scale implies 
for the quark-gluon correlations probed by the twist-4 
and 3 operators. 

A microscopic model incorporating the short-distance 
scale associated with chiral symmetry breaking is the in- 
stanton vacuum, in which the quarks obtain a dynamical 
mass, M, by interaction with a "medium" of instantons 
of characteristic size n k, 0.3 fm, see Refs. Hi m for 
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a review. The fundamental assumption of diluteness of 
the instanton medium (smallness of the packing fraction) 
imphes that, parametrically, M <C p^^ ■ Since the range 
of the chiral forces binding the constituent quarks is of 
the order M~^, this hierarchy translates into a "two- 
scale picture" of hadron structure, Eq. Q. The higher- 
twist matrix elements Q and (jSJ were calculated in this 
model in Refs. 17] . It was found that the flavor- 
nonsinglet twist-4 matrix element is determined by the 
inverse instanton size {i.e., the inverse size of the con- 
stituent quark), f2~'^ and has a sizable nega- 
tive value. The twist-3 matrix element, d2, however, 
is parametrically suppressed, d2 ^ (Mp)^, in agreement 
with the experimental data (this prediction was later con- 
firmed by lattice QCD calculations |0|). An interesting 
question is whether these results depend specifically on 
the assumption of instantons as the dominant gluonic 
vacuum fluctuations, or whether they already follow from 
the more general "two-scale picture" of hadron structure. 



In this paper we further explore the connection be- 
tween the short-distance scale due to chiral symmetry 
breaking and the twist-4 and 3 quark-gluon correlations 
governing the corrections to the nucleon spin struc- 
ture functions. In Sec. ^ we argue that, on general 
grounds, the flavor-nonsinglet twist-4 matrix element, 
f2 '^, has a sizable negative value of the order p^^, due 
to the presence of sea quarks with virtualities ~ p~^ in 
the nucleon wave function. The twist-3 matrix element, 
(i2, however, is not related to the scale p~^ . Our argu- 
ments provide additional insights into the origin of the 
instanton vacuum results, and suggest that they may be 
of more general nature. In Sec. IIIII we present the re- 
sults of a next-to-leading order (NLO) QCD analysis 
of the world data for g\ and qV" , including recent data 
from the Jefferson Lab HaU A and COMPASS 
experiments, in which we extract the flavor-nonsinglet 



twist-4 matrix element, /^~''. The sign and order-of- 
magnitude are found to be in agreement with the our 
qualitative arguments (as well as the instanton vacuum 
results). This analysis extends previous phenomenolog- 
ical estimates of the higher-twist contribution to polar- 
ized deep-inelastic scattering by Leader et al. 0, 
In contrast to other analyses [ij we perform the QCD 
fit to the data (including higher-twist corrections) for 
the x-dependent structure function, computing the mo- 
ments only at the last stage, by integration of the fit. 
Furthermore, we extract directly the isovector (flavor- 
nonsinglet) twist-4 matrix element, which is practically 
scheme-independent and thus provides a much cleaner 
probe of the non-perturbative quark-gluon correlations 
than the singlet matrix element. Finally, in Sec. IIVI we 
comment on the implications of the "two-scale picture" 
of hadron structure for the x-dependence of the higher- 
twist contribution and quark-hadron duality in the spin 
structure functions. 



II. HIGHER-TWIST MATRIX ELEMENTS AND 
THE SIZE OF THE CONSTITUENT QUARK 

The spontaneous breaking of chiral symmetry implies 
the existence of a short-distance scale in the QCD vac- 
uum, significantly smaller than the typical hadronic size. 
It is determined by the characteristic size, p, of the non- 
perturbative gluon field configurations which change the 
quark chirality, see Fig. ^ A gauge-invariant measure of 
this size is the ratio of the dimension~5 "mixed" quark- 
gluon condensate to the usual dimension-3 quark con- 
densate, 



(4) 



Here cr°^ = (i/2)[7", 7^], and Va = 9a - igA^ is the 
covariant derivative. The two representations of the 
dimension-5 operator are related by the identity 



(5) 



and the Heisenberg equations of motion for the quark 
fields. Lattice simulations indicate that rriq > 1 GeV^ at 
a normalization point of /i 1 GeV 0, J2a| (substan- 
tially larger values were obtained in Ref. '23|). A more 
precise interpretation of the ratio Q in terms of a size 
of field configurations becomes possible with specific as- 
sumptions about the shape of these configurations in a 
given gauge. In the instanton vacuum where the 

chirality-flipping field configurations are (anti-) instan- 
tons in singular gauge, one has 



Ap- 



(6) 



at the scale (jl '-^ [2J|. The lattice results are con- 
sistent with an average instanton size of p « 0.3 fm. We 
note that in the presence of more than one light quark 
fiavor, the chirality flip due to the instanton happens in 
a many-fermionic interaction ('t Hooft vertex), whence 
the scale p may also be interpreted as the range of chiral- 
symmetry breaking quark-quark interactions in the QCD 
vacuum. 

A simple heuristic argument suggests that the twist-4 
matrix element, /2, Eq. may be related to the short- 
distance scale, p. Substituting in the twist-4 operator the 




FIG. 1; Origin of the short-distance scale, p, associated with 
cliiral symmetry breaking. 
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FIG. 3: The "sea" quark contribution to the axial charge, 
giving rise to the large-virtuality tail {—k^ ~ P~^) in the 
distribution f{k), Eqs. ^ and (IHl . 



FIG. 2: Distribution of quark virtuaUties, , in the proton's 
isovector axial charge, as given by the effective chiral theory, 
Eg. llH . The twist-4 matrix element, /^~'*, is of the order of 
the integral of k^ f{k). 



gluon field by the commutator of covariant derivatives, 
Eq. (jSJ , and making use of gamma matrix identities and 
the Heisenberg equations of motion of the quark fields, 
one can convert the twist-4 operator to the form 



/2 : V57''75(-V2)V. 



(7) 



In this form, it can be compared with the axial current 
operator, which measures the quark contribution to the 
nucleon spin. 
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(8) 



We see that the operator ^ measures the correlation of 
the spin of the quarks with the square of their canon- 
ical momentum. The existence of the short-distance 
scale, p, implies that, generally speaking, nucleon cor- 
relation functions involve "sea" quarks with virtualities 
(four-momenta squared) up to the scale . If the spin 
of these quarks is correlated with the nucleon spin, one 
would expect the flavor-nonsinglet twist-4 matrix ele- 
ment to be of the order 



/2 
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(9) 



Here we limit ourselves to the flavor-nonsinglet (u — d) 
operators, which are not affected by the U(1)a anomaly 
(henceforth, gA denotes the nucleon's isovector axial cou- 
pling). Furthermore, if the flavor-nonsinglet sea quarks 
were polarized along the direction of the nucleon spin, 
and if the quark canonical momentum squared on aver- 
age takes positive values, as suggested by Eq. one 
would conclude that /^"'' < 0. 

The above argument supposes the existence of a dy- 
namical mechanism which polarizes the flavor-nonsinglet 
"sea" quarks in the nucleon. In fact, such a mechanism 
can be found in the effective dynamics resulting from the 
spontaneous breaking of chiral symmetry. In the large- 
Nc limit of QCD, the effective dynamics at distances ^ p 
can be described by massive "constituent" quarks, cou- 
pled to a Goldstone pion field in a chirally invariant way 



^cff 



(10) 



The dynamical quark mass, and the coupling to the pion 
field, are active for quark virtualities < p~^, whence 
p ^ can be interpreted as the "size" of the con- 

stituent quark. The large- A^c limit implies a semiclassical 
description of the nucleon, in which the nucleon is char- 
acterized by a classical pion field of size R ~ |27j . 
Thus, the effective dynamics described by Eq. ifTUIl em- 
bodies the "two-scale picture" of the nucleon described 
in the introduction, Eq. 

Consider now the matrix element of the flavor-non- 
singlet axial current operator, i.e., the nucleon's isovec- 
tor axial charge, in the effective theory of constituent 
quarks coupled to a pion field, Eq. (^UJ). After ap- 
propriate projections on nucleon states with definite 
spin/flavor/momcntum quantum numbers, the matrix el- 
ement can be represented in the form 
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i2n) 



(11) 



where k represents the four-momentum of the quark con- 
nected to the axial current operator. We are interested in 
the distribution of quark virtualities, fc^, in this integral 
|50| . This distribution consists of two components, shown 
schematically in Fig. |21 The bulk of the nucleon's axial 
charge is carried by valence and sea quarks with virtual- 
ities of the order of the size of the nucleon, — fc^ ^ R^^. 
In addition, however, there is a contribution from "sea" 
quarks interacting with the classical pion field, which 
have virtualities extending up to —k^ ~ P~^- Assum- 
ing that p R, the power behavior and the coefficient 
of the this large-virtuality "tail" are completely deter- 
mined by general features of the chiral dynamics. They 
follow from the leading term in the long-wavelength ex- 
pansion (gradient expansion) of the quark loop in the 
background pion field (see Fig. |3Jl : 



/(fc) 
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(fc2)2 

where the coefficient is given by 



(12) 
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d'^x tr 



iT''Ul{x)daU,i{x) . (13) 



Here Uc\{x) = exp [z7r"[(a;)T"] denotes the static classical 
pion field in the nucleon rest frame. In fact, the ex- 
pression Hi;-{|l coincides with the leading-order gradient 
expansion of the axial charge induced by the classical 
pion field (as used e.g. for calculating gA in the Skyrme 
model). The power behavior H12(l implies that the in- 
tegral depends only logarithmically on the short- 
distance scale p: 

gA ~ log-^. (14) 

Thus, gA is dominated by quark virtualities of the order 
<C p~^, and the presence of the large-virtuality tail 
in the distribution of Fig.|21is of minor importance. 

The presence of the large-virtuality tail in the ax- 
ial charge distribution becomes crucial, however, for the 
twist-4 matrix element, f2~'^- When passing from QCD 
to the effective theory of constituent quarks, QCD oper- 
ators must be "translated" to operators in the effective 
theory. It is natural to assume that in the case of the 
twist-4 QCD operator Q the translation is given by 

^7^75(-V2)V|qcd C^7^75(-52)V'|cff, (15) 

where C > is a coefficient of order unity. This was 
explicitly demonstrated in the instanton vacuum model, 
where the operator matching follows from the integra- 
tion over instanton-type gauge field configurations |l6|. 
By analogy with the axial current, the matrix element 
of the constituent quark operator 1)15(1 is now given by 
the integral of times the momentum distribution of 
gA, Eq. (fTT|l ). This factor suppresses contributions from 
virtualities —k^ ~ R~^, and enhances the contributions 
from the large-virtuality tail. As a result, the twist-4 
matrix element is parametrically of the order 

fr" - p-\ (16) 

with the coefficient proportional to gA, in agreement with 
our estimate ©■ Furthermore, since the quark virtuali- 
ties in the "tail" of the distribution are spacelikc, k^ < 0, 
we conclude that /;^~'' < 0. 

The integral determining f2 '^ in the effective chiral 
theory contains a would-be quadratic divergence, which 
is cut off at the scale p~^. While the parametric order 
and sign of f2~'^ follow from general features of the chiral 
dynamics, the numerical value is very sensitive to the 
precise way in which the UV cutoff (i.e., the finite size 
of the constituent quark) is implemented. The instanton 
vacuum model, which implies a definite form of the UV 
cutoff, resulting from the fermionic zero modes of the 
instantons |26l |. and which also allows one to uniquely 
determine the matching coefficient for the effective twist- 
4 operator, Eq. (jTs)) [ig], gives (with p = 0.3 fm) Tt^ 

f^-'^ « -O.bgAP'^ = -0.22 GeV^ (17) 

However, in view of the principal uncertainties in the 
modeling of the dynamics of quark field modes of vir- 
tualities p^^ this result should be viewed as a rough 
estimate. 
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f^-'' [GeV^] 


Instantons [16, 17] 


-0.22 


QCD sum rules (Balitsky et al.) [2£] 


-0.16 ±0.04 


QCD sum rules (Stein et al.) [301 


-0.06 ±0.02 


Bag model \5j 


±0.09 



TABLE I: Comparison of theoretical estimates of the flavor- 
nonsinglet twist-4 proton matrix element, f2~'^- 



On general grounds, one expects that in QCD with 
cutoff regularization the matrix elements of twist-4 op- 
erators are proportional to the square of the cutoff |28j| : 
this is realized in our approach if the cutoff is identi- 
fied with the scale p^^. We do not consider here the 
logarithmic scale dependence of the twist-4 matrix ele- 
ment which results from the standard composite operator 
renormalization Q ; this dependence is much weaker than 
the principal uncertainty in our estimates of the twist-4 
matrix element. 

It is interesting that the estimates of f2~'^ obtained 
in QCD sum rule calculations |2^|33| agree in sign and 
order-of-magnitude with our estimate, and with the in- 
stanton vacuum result, see Tabled (for a critical discus- 
sion of these calculations, see Ref. j3l|). Our estimate 
disagrees with the bag model 0, which gives a positive 
result for /^^'^ (note the different convention for the sign 
of /2 in that paper). This model, however, does not re- 
spect the QCD equations of motion, and therefore cannot 
claim to give a realistic description of quark-gluon cor- 
relations in the nucleon. 

When applying the same reasoning as above to the 
twist-3 matrix element, ^2, we find that after the sub- 
stitution (O the quark-gluon operator does not produce 
a contracted covariant derivative. In this operator, all 
derivatives are "kinematic", i.e., they are needed to sup- 
port the spin of the matrix element. This operator does 
not probe the virtuality of the quarks in the nucleon, 
and its matrix element does not receive essential contri- 
butions from quark virtualities of the order p~^. Thus, 
the correlations probed by the twist-3 matrix element 
are of essentially different nature as those probed by the 
twist-4 one. Beyond this qualitative difference, we see no 
simple way to estimate the twist-3 matrix element more 
accurately on grounds of general features of the effective 
dynamics alone (see also Ref. |32|). In the instanton vac- 
uum, an additional suppression of d2 results from the fact 
that the coefficient of the corresponding effective opera- 
tor in the effective chiral theory is parametrically small 
in the instanton packing fraction (i.e., d2 is zero in the 
single instanton approximation). This suppression ap- 
pears to be due to the 0(4) invariance (in the Euclidean 
metric) of the instanton field |16|. In a sense, the twist-3 
operator is a much more subtle probe of non-perturbative 
quark-gluon correlations in the QCD vacuum than the 
twist-4 operator, whose matrix element can be estimated 
on general grounds. 
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III. TWIST-4 MATRIX ELEMENT FROM A 
QCD FIT TO POLARIZED DIS DATA 

It is interesting to see to which extent our quahtative 
conclusions about the higher-twist matrix elements are 
supported by the experimental data. To this end, we 
attempt to extract the flavor-nonsinglet twist-4 matrix 
element, /^~'^, from the power corrections to the spin 
structure functions and g". To accuracy the 
tree-level QCD expansion for the first moment of — 
is given by [1110 
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A/2 

12 "2 



27Q2 



u — d 



+fr') (18) 



(when QCD radiative corrections are included, the coef- 
ficients acquire a logarithmic Q^-dependence) . The first 
term is the leading-twist (LT) contribution, proportional 
to the matrix element of the flavor-nonsinglet twist-2 
spin-f operator, Oq"'', with Oq"'' = gA (Bjorken sum 
rule). The second term represents the target mass correc- 
tions (TMC), proportional to the spin-3 twist-2 matrix 
element, Oj"''. The third term is the dynamical higher 
twist (HT) contribution, involving the twist-3 and 4 ma- 
trix elements Q and (Q. The twist-3 matrix element, 
d2, has been extracted from independent measurements 
of the third moment of the spin structure function 52 
(with the Wandzura-Wilczek contribution subtracted). 
The SLAC E155X experiment and the recent Jeffer- 
son Lab Hall A analysis W\ report values of " ^ 10"^, 
in g ood agreement with the instanton vacuum prediction 
[T6| . With these values, the contribution of the Mj^d'^^'^ 
term to the corrections in (|18|) is more than an 

order of magnitude smaller than the instanton estimate 
for f2~'^, Eq. (|17(l . This supports the qualitative con- 
clusion from our two-scale picture, that the dominant 
power corrections are those associated with the short- 
distance scale, p. We shall thus neglect the M]qd^~'^ 
term compared to f2~'^ in Eq. H18|) . and ascribe the phe- 
nomenological power correction entirely to the twist-4 
matrix element. This theoretical simplification will be 
justified a posteriori by the fact that the numerical value 
of f2 '^ extracted in this way is indeed much larger than 

the measured M^dj"*^- 

The dynamical higher-twist contribution to the x~ 
dependent structure functions, 5i(a;,Q^) and gi{x,Q'^), 
has been extracted from NLO QCD fits to the world data 
(see references in Ref. 33]), including the new data 
from the Jefferson Lab Hall A experiment [T^ , as well as 
the deuteron data from COMPASS 113. These fits are 
based on the ansatz 

9rM^) = 5r(^,Q')LT + TMC + (19) 

where the leading-twist contribution (including target 
mass corrections) is given by the Leader-Stamenov- 
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FIG. 4: The difference of higher twist corrections to the 
proton and neutron spin structure functions, h^{x) — h"{x), 
Eq. 119II . as extracted from our NLO QCD fit to the world 
data (see references in Ref. [s^), including recent (/" data 
from the Jefferson Lab Hall A experiment [l9|. as well as 
deuteron data from COMPASS IIS. 



Sidorov parametrization of the polarized parton densi- 
ties [s^, and parametrizes the dynamical higher- 
twist corrections. In order to extract directly the flavor- 
nonsinglet higher-twist correction, we have modified 
the fit procedure of Refs. H [s^ and parametrized not 
}iPi"-)(^x) individually, but their difference and sum. The 
kinematic cuts applied in the new fit are the same as in 
the fit of Ref. 0, > 1 GeV^ and W^>A GeV^. The 
results for the difference, h'P{x) — h'^{x), obtained in this 
way is shown in Fig. 0] 

Integrating the higher-twist contribution over the x— 
range covered in our fit we get 



0.75 



dx [hP{x)-h''{x)] 



-0.046 ± 0.016 GeV^ (20) 



0.017 



If we neglected the contribution from the large-a; region 
(we shall argue in Section Hvl that this is consistent with 
our two-scale picture of the structure of the nucleon) , as 
well as the small-x region, and regarded the integral H20() 
as an estimate of the first moment, we would obtain [cf. 
Eq. IIHI)] 



rv 
J2 



-0.31 ±0.11 GeV" 



(21) 



For a rough data-based estimate of the contribution of 
the large-x region, we consider the integrals up to x = 1 
computed with a linear extrapolation of the fit into this 
region, based on the two rightmost data points in Fig. ^ 
and a constant extrapolation, based on the rightmost 
data point only (we neglect the small-x contribution). 
Taking the average of the two integrals as an estimate of 
the central value, and half the difference as an estimate 
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of the error, we get 

[ dx [hP{x) - h"{x)] = -0.028 ± 0.019 GeV^, (22) 
Jo 

corresponding to 

f^-'^ = -0.20 ± 0.14 GeV^ (23) 

Both estimates (|21f) and (|23|l agree in sign and in order- 
of-magnitude with our quahtative prediction, Eq. 
and the resuh obtained in the instanton vacuum, 
Eq. |T7|). 

In the QCD fit, the leading-twist parton densities and 
coefficient functions were taken in the AIS scheme. The 
studies of Ref. ^ found that the higher-twist corrections 
to gi in different factorization schemes (in particular, in 
the MS and JET schemes) coincide within errors. Note 
also that Eq. H21|) agrees well with our previous estimate 
based on the JET scheme [s^. In that estimate, also pre- 
liminary HERMES data on the deuteron structure func- 
tion were taken into account The recent analysis of 
Ref. shows that inclusion of the new HERMES pro- 
ton and deuteron data 0| does not substantially change 
the higher-twist contribution, whence we have not ex- 
plicitly included these data in our present analysis. 

Our result for /^^"'^ agrees well with that obtained by 
Deur et al. [T]| in a recent analysis of power corrections 
to the Bjorken sum rule (their = ^f2~'^ in our 

conventions). Our method differs from that of Ref. 0| 
in that we perform the QCD fit to the data (including 
higher-twist corrections) for the x-dependent structure 
function, computing the moments only at the last stage, 
by integration of the fit. Nevertheless, the results for the 
higher-twist contribution from both approaches are com- 
parable, which is very encouraging. Our result disagrees 
in sign with that obtained in an earlier analysis of the 
Q^-dependence of the Bjorken sum rule (33|, which com- 
bined an empirical parametrization of the electroproduc- 
tion cross section in the resonance region with the QCD 
parametrization in the DIS region, see also Refs. |3i\i40h |. 

IV. DISCUSSION AND OUTLOOK 

In our discussions so far we considered the twist-4 cor- 
rection to the first moment of gi{x), which is related to 
the matrix element of the local twist-4, spin-1 opera- 
tor, /2. Much more information is contained in the x- 
dependence of the higher-twist contribution. A detailed 
study of the x-dependence of the twist-4 contribution on 
the basis of chiral dynamics and the instanton vacuum is 
beyond the scope of the present paper. Here we would 
like to offer only some general comments on this problem. 

The idea of a two-scale picture of the nucleon out- 
lined in Section ^ cf. Eq. Q can also be expressed in 
a partonic language. In this formulation, constituent 
quarks/ antiquarks appear as correlations in the trans- 
verse spatial distribution of quarks/ antiquarks and glu- 
ons, with a transverse size, p, significantly smaller than 



the transverse size of the fast-moving nucleon, R. (This 
formulation is in fact used to discuss the effect of the con- 
stituent quark structure of the nucleon on high-energy 
pp scattering with multiple hard processes 41J.) This 
picture applies to average values of the quark/antiquark 
longitudinal momentum fraction, x < 0.5, for which the 
transverse size of configurations in the nucleon wave func- 
tion is of the order of the typical hadronic size. In this 
formulation, the twist-4 correction to the spin struc- 
ture function can be related to the average transverse 
momentum squared of the polarized quarks/antiquarks 
in the nucleon l43l |. Extending the reasoning of 
Section ^1 one would argue that due to chiral dynam- 
ics the flavor-nonsinglet polarized "sea" quark distribu- 
tion involves transverse momenta squared of the order 
/c^ ~ p~^. This suggests an interesting connection be- 
tween the flavor-nonsinglet twist-4 corrections to 171 and 
the large flavor asymmetry of the twist-2 sea quark dis- 
tribution, Au{x) — Ad(x) predicted by chiral dynamics 
in the large- A'^c limit |4J|. The result of our QCD fit. 
Fig. ^ indicates that the twist-4 correction is indeed lo- 
calized at relatively small values of a; ~ 0.1, supporting 
the connection with the sea quark distribution. A similar 
connection between the flavor asymmetry and the pres- 
ence of large transverse momenta in the sea quark distri- 
butions was noted in Ref. "45], where instanton-induced 
sea quark components in the nucleon wave function were 
considered in a phenomenological model with no refer- 
ence to the large- Ac limit. 

The existence of the short-distance scale due to chi- 
ral symmetry breaking also has some interesting qualita- 
tive implications for quark-hadron duality in polarized 
DIS. The two-scale picture of the structure of the nu- 
cleon, cf. Eq. ||2Jl, implies a parametric classification of 
the hadronic excitation spectrum of the nucleon. Nu- 
cleon resonances such as the A are excitations with en- 
ergies (in the CM frame) of the order E ^ R^^. They 
correspond to changes of the state of motion of the con- 
stituent quarks over distances R, which do not affect 
the internal structure of the constituent quark at dis- 
tances ^ p. Excitations of energy E ^ p^"^ belong to 
the non-resonant hadronic continuum. Switching to the 
quark language, our arguments of Section ^ show that 
the twist-2 quark distribution (here, the axial coupling, 
qa) arises mainly from field configurations with ener- 
gies/momenta of the order R^^ , while the twist-4 quark- 
gluon correlations are dominated by energies/momenta of 
the order p~^ . Comparing the hadronic and the quark de- 
scription, we conclude that quark-hadron duality should 
"work" for the twist-2 part of the structure function (i.e., 
the Q^-independent part) already when summing over 
hadronic excitations with energies E ~ but for the 
twist-2 plus twist~4 part of the structure function (i.e., 
to accuracy 1/Q^) only when summing over hadronic ex- 
citations with energies E ^ p~^. In practice, this means 
that quark-hadron duality in the structure function to 
accuracy may require integration over a signifi- 

cantly larger duality interval than duality to accuracy 
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(Q^)". This needs to be taken into account when trying 
to extract higher-twist matrix elements from resonance- 
based parametrizations of the structure functions. To 
summarize, the two-scale picture makes a clear para- 
metric distinction between resonance and higher-twist 
contributions to the structure function. This qualitative 
prediction is supported by the fact that the phenomeno- 
logical twist-4 contribution to gi (see Fig. 0J| seems to 
be dominated by small values of a;, below the resonance 
region. 

The two-scale picture described here is close in spirit to 
the loffe-Burkert parametrization of the dependence 
of the first moment of gi , in which the contribution 
from the A resonance is separated from the continuum, 
and the leading power corrections are associated with 
the continuum contribution. The characteristic mass 
scale governing the power corrections in this parametriza- 
tion, /i^ = Mp, is numerically close to value associated 
with the constituent quark size, p^^ ^ (0.3fm)~^ = 
(600 MeV)^. Also, the analogous parametrization for the 
Bjorken sum rule (p — n) gives negative sign of the twist- 
4 correction, in agreement with our qualitative predic- 
tion. However, the numerical value of the twist-4 correc- 
tion obtained from the loffe-Burkert parametrization is 
substantially larger than the instanton vacuum estimate, 
corresponding to f^'"' « -2.3 gA A*^ = 1-7 GeV^. 

More generally, the two-scale picture of hadron struc- 
ture allows one to draw some conclusions about global 
properties of the transition from high to low in the 
nucleon spin structure functions (i.e., going beyond the 
leading corrections). Since the characteristic mass 

scale for the power corrections is set by the size of the 
constituent quark, one should expect the twist expansion 
to break down at momenta of the order ^ p^^ . For 
the extraction of the leading (1/Q^-) corrections from 
QCD fits to the data this implies that one should restrict 
oneself to the range ^ where the leading term in 
the series dominates (in our fit presented in Section IlIII 



g2 > IGeV^). 

The arguments presented in this paper can also be ex- 
tended to higher-twist corrections to unpolarized deep- 
inelastic scattering. A new feature in the unpolarized 
case is the appearance of twist-4 operators measuring 
quark-quark correlations ("four-quark operators"), in 
addition to quark-gluon correlations of the type encoun- 
tered in the polarized structure functions. Within our 
two-scale picture, the quark-quark correlation matrix el- 
ements are of the order R~^, and thus paramctrically 
suppressed compared to the quark-gluon ones, which are 
of the order p^^ . This qualitative conclusion seems to 
be in agreement with a joint analysis of the twist-4 cor- 
rections to F2 and Fl, see Ref. (48^ for details. The role 
of the size of the constituent quark in power corrections 
to unpolarized structure function moments was also dis- 
cussed in a different approach in Ref. 47] . 

To summarize, we have argued that the leading power 
corrections to the nucleon spin structure functions are 
governed by the short-distance scale due to the spon- 
taneous breaking of chiral symmetry — the size of the 
constituent quark. The qualitative statements following 
from this assumption are supported by the result of a 
QCD fit to the present polarized DIS data. The argu- 
ments presented here may eventually serve as the basis 
for an "interpolating" model of the nucleon spin struc- 
ture functions, connecting the scaling region at large 
with the photoproduction point. 
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the coefficient of the twist-4 matrix element, /2, Eq. Q, 
such that it has dimension (mass)^, i.e., we do not make 
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mass). We shall see below that the value of /2 is deter- 
mined by a physical mass scale unrelated to the nucleon 
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